Transition from Bose-Einstein condensate to Berezinskii-Kosterlitz-Thouless phase 
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We obtain a phase diagram for a trapped two-dimensional ultra-cold Bose gas. We find a critical 
temperature above which the free energy of a state with a pair of vortices of opposite circulation is 
lower than the one for a vortex-free Bose-Einstein condensed ground state. We identify three distinct 
phases which are, in order of increasing temperature; a phase coherent Bose-Einstein condensate, a 
vortex pair plasma with a fluctuating condensate phase, and a thermal Bose gas. The existence of 
the vortex pair phase could be verified using current experimental setups. 
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Introduction 

The theory of phase transitions is fundamental to the 
statistical mechanical description of nature P|. Techni- 
cally, phase transitions involve non-analyticity in the free 
energy with a thermodynamic variable. The solid-liquid- 
gas transitions in water are common examples of first- 
order phase transitions, involving latent heat. In second- 
order phase transitions, such as the ferromagnetic transi- 
tion, latent heat is absent at the critical point. If, in ad- 
dition, the transition does not involve a symmetry break- 
ing, it is classified as of infinite order. A famous example 
of such is the Berezinskii-Kosterlitz-Thouless (BKT) vor- 
tex unbindingtransition in two-dimensional superfiuid 
systems |3, IMQ ■ Evidence for this phase transition has 
been experimentally demonstrated in liquid helium thin 
films j^, in superconducting Josephson-junction arrays 
0, and in spin-polarized atomic hydrogen 0. 

Bose-Einstein condensation (BEC) is a purely 
quantum-statistical second-order phase transition occur- 
ring even in the absence of interactions 0, 0, Since 
the realization of BEC in 3D harmonic traps [ill ll^ . 
there has been an increased interest and debate on the 
nature of the normal-superfluid phase transition in the 
2D Bose gas. It has long been established that BEC 
does not occur in homogeneous two-dimensional systems 
in the thermodynamic limit. This is formally proved in 
the Mermin-Wagner-Hohenberg theorem |l3l [iJI . Nev- 
ertheless, an external confining potential modifies the 
density of states in such a manner that, in a 2D ideal 
Bose gas, condensation into a single quantum state is 
predicted to occur at finite temperature fl^ . It has been 
shown, however, that long-wavelength phase fluctuations 
may destroy the global phase coherence of the condensate 
below the ideal gas critical temperature, Tq, for BEC. 
The condensate can then only exist in the form of a so- 
called quasi-condensate E El El El HO, El [13 ■ The 
connection between such a quasi-condensate and a BKT 
vortex-pair plasma remains unclear. Suggested methods 
for achieving the required two-dimensionality for observ- 
ing the BKT phase transition in dilute atomic gases in- 
clude application of tight axial trapping E^l , use of op- 
tical lattice potentials 13, and rapidly rotating gases 




FIG. 1: Computed spatial distribution of the squared mod- 
ulus of the order parameter containing a pair of vortices of 
opposite circulation. The arrows represent the local super- 
fiuid currents. The vortex pair shown is widely separated for 
demonstrational purposes. 



with hard- wall radial confinement [25j . In this Letter we 
show, based on a free energy argument, that a 2D har- 
monically trapped Bose gas may exhibit both the BEC 
and BKT-like superfiuid phases, separated by what we 
believe to be a first-order phase transition. At higher 
temperature, the latter is expected to undergo a vortex 
unbinding transition to a normal state analogous to the 
usual BKT transition for the homogeneous system. 

Free energy 

We consider a dilute gas consisting of N bosonic par- 
ticles of mass m spatially confined in two dimensions by 



an external harmonic potential ^(r) 



■/2, where 



ij^ is the radial trap frequency. We assume the existence 
of a macroscopic ground state condensate wavefunction 
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■0(r), which also serves to describe the order parameter 
of the system in the BEC phase. Our aim is to calculate 
and compare the Hclmholtz free energy F ^ E ~ TS 
for a phase coherent ground state and for an excited 
state containing a vortex pair in the ordering field, see 
Fig.n for various configurations. Due to entropy consid- 
erations, one expects to find a critical temperature, Tc, 
above which a thermally activated transition to the state 
containing vortex-pair excitations becomes favourable. 

Isolating the contribution due to the condensate, 
Eq{T), the total internal energy of the system may be 
written as E{T) = Eo{T) + E{T). The required conden- 
sate energy Eq{T) is determined by the functional 

E,{T)^ J (^|l|VV'(r)p + y(r)|V'(r)p + ||V'(r)|4^ dr 

(1) 

where g is the constant coupling parameter for the parti- 
cle interactions, and Nq = / |-0(r)pdr denotes the num- 
ber of particles in the condensate. In order to implicitly 
incorporate the temperature dependence in Eo{T), we 
employ the 2D ideal-gas result for the condensate frac- 
tion 



N \To 



(2) 



Alternatively, the same where fc^To = huj i_^J N / C,{2) and 
C is the Riemann zeta function. The shift in the critical 
temperature from the 2D ideal-gas prediction Tq is only of 
a few percent, due to the weak particle interactions |22l |. 
We therefore expect the above approximation to yield a 
good estimate for the pair creation transition tempera- 
ture in the system considered here. The internal energies 
for topologically distinct order parameter configurations 
are calculated both analytically within the hydrodynamic 
approximation and computationally as described along 
with the numerical methods. 

To evaluate the free energy difference, one also needs 
to calculate the entropy change A5' = fc^lniy, due 
to the multiplicity of order parameter configurations 
containing a pair of vortices. The statistical weight, 
W = , is obtained by allowing one vortex to 

reside anywhere within the Thomas-Fermi radius Rtf ~ 
yj2ii/mijj'\_ of the condensate, where fji is its chemical po- 
tential. The partner vortex then has 27r available nearest 
neighbour sites, assuming that the pair is closely bound 
when it forms. Approximating the radius of the area 
occupied by a vortex to be equal to the healing length, 
^ = /2m^, of the order parameter, we obtain the 
configurational entropy AS{T) = ks In (87r/x^/?i^w^) due 
to the available microstates containing a pair of vortices. 
Here we assume that all of these microstates possess the 
same internal energy although in a strict sense this is 
not true, since the gas is inhomogeneous. However, mi- 
crostates of equal energy with the vortex pair residing 
outside the condensate centre may be chosen to be ones 



with correspondingly larger vortex pair separation, bal- 
ancing the energy reduction due to the inhomogeneous 
superfluid density. By applying the Thomas-Fermi ap- 
proximation for the chemical potential, the relevant free 
energy difference per particle is given by the formula 



AF(T) AEo{T)No T In {8 gNom/h^) 



N 



No N To VC(2)iV 



(3) 



Guided by the results for the excitation energy of a single 
vortex at finite temperatures |2^ , we have assumed that 
the energy difference associated with uncondensed atoms 
AE{T) < AEoiT) in Eq. ®. In the same fashion, the 
entropy difference due to the quasiparticlc excitations is 
neglected. 

The critical temperature, Tc, for the pair creation tran- 
sition is determined by the condition AF{Tc) = 0. Let- 
ting r+ and r~ denote the positions of the oppositely 
circulating vortices, we define the energy of a pair exci- 
tation per condensate particle a(r+,r~) = AEo{T)/No- 
Hence, the transition temperature is implicitly given by 
the formula 



Tc 
To 



No 



huj± y/Nln {8gNom/h?)' 



(4) 



In the hydrodynamic approximation the energy of a 
vortex pair may be estimated to be [l(l | 



"—\n- + 2Ec, 

m 4 



(5) 



where no is the superfluid particle density, D — 2S, is the 
pair separation, and E^ denotes the energy associated 
with a single vortex core. Assuming that the condensate 
density within an isolated vortex core increases quadrat- 
ically as a function of the distance from its centre, we 
obtain 



E, 



vp 



2-K^h^ 

gm 



In 2 



(6) 



for a vortex pair located in the vicinity of the trap centre. 
Inserting Eq. © for AEo{T) in a, one finds that the 
critical temperature is given by 



To 



P + g\n{8gmNo/h? 



(7) 



where, (3 = 2T:'^h^{\n2 + l/2f/im. In the limit of the 
non-interacting ideal-gas limg_,o(T'c/To) = 1 and the 
BKT phase vanishes. In homogeneous systems the BEC 
phase does not exist — a result that can be recovered by 
considering the large g limit in which case the condensate 
density becomes homogeneous on the length scale of the 
vortex core. This conclusion is also reached if the ther- 
modynamic limit, w — > and N ^ oo with Nuj'^ held 
constant, is taken [loj . 
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FIG. 2: Solution of Eq. llOt as a function of the total particle 
number A*". The phase boundaries are found by applying dif- 
ferent estimates for the vortex pair energy a. The diamonds 
correspond to the full numerical calculation whereas the con- 
tinuous lines are plotted using Eq. © with (solid line), or 
without (dashed line) the core energy contribution. 



Numerical results 

In order to assess the validity of the hydrodynamic 
approximation we calculate the pair excitation energy a 
by solving numerically the Gross-Pitaevskii equation 



2m 



V' + V{r)+g\^{r)\']^{r)^fiiP{r) (8) 



for a variety of order parameter configurations. The vor- 
tex pair states are generated by initially phase imprint- 
ing the two vortices with opposite circulation in the or- 
der parameter field as depicted in Fig. ^ The subse- 
quent damped real-time propagation of V'(r) under the 
Hamiltonian corresponding to the Ihs of Eq. © serves to 
provide an ensemble of variational wavefunctions whose 
energies are given by the energy functional, Eq. 

Computationally, all dimensional parameters in Eq. (jS)) 
can be embedded into a single constant C = gNQ/hio±aQ 
where qq = y^h/2muj± defines the characteristic length 
scale in the system. We model the experimentally achiev- 
able quasi-2D systems and hence the coupling constant 
g = y/Snh'^a/maz, where a is the s-wave scattering 
length, Qz = y^h/mujz, and uJz is the axial trapping 
frequency. We have also verified that using this quasi- 
2D value for the particle interactions is indistinguishable 
from the full T-matrix result [itI Erj . The equation for 
C may thus be expressed as 



C = 4V2^N- 

a-z \ \Tq, 



(9) 



explore the full parameter space of the problem. Thus, 
once the vortex pair excitation energy a is calculated, 
the critical temperature is obtained from the equation 



Tc 
To 



Nn 



In {16V2^Noa/az)' 



(10) 



Modelling the 2D experimental system at Oxford |23l |. 
we choose {co±,ujz} ^ 2tt x {3,2200} Hz. We have cal- 
culated the pair excitation energy a for a range of pa- 
rameters C, r+ and r~. This energy depends on the 
separation |r+ — r^| of the vortices due to their logarith- 
mic attractive interaction and is proportional to the local 
condensate density at the vortex locations. The main dif- 
ficulty in obtaining accurate estimates for the energies of 
the vortex-pair states is the fact that well defined vortex- 
pair eigenstates do not exist. Additional problems arise 
with large condensates for which the vortex cores are 
very small and the energy of the ground state is many 
orders of magnitude larger than that of the vortex pair 
excitation. 

The solution of Eq. (fTn|l is plotted in Fig. ISJ as a 
function of the total particle number N, using the Ox- 
ford trap parameters and three different estimates for the 
vortex pair energy a. The diamonds correspond to the 
full numerical calculation whereas the continuous lines 
are plotted using Eq. © with (sohd fine), or without 
(dashed line) the core energy contribution. 

The free energy difference changes sign at the phase 
boundary and vortex pair creation becomes favourable at 
higher temperatures. It is to be noted that in a real sys- 
tem the thermally activated vortex pair nucleation would 
be likely to occur first near the low density perimeter of 
the condensate and in this respect the estimated phase 
boundaries in Fig. Q correspond to upper bounds on 
the true critical temperature. 
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and therefore, by altering the single parameter C, one can 



FIG. 3; Phase diagram for a 2D trapped Bose gas at ultra-low 
temperatures. 



4 



Thus we arrive at the phase diagram shown in Fig. (O 
for the dilute 2D gas of bosonic particles at ultra-low 
temperatures. For the lowest temperatures the gas ex- 
ists in a Bose-Einstein condensed state. With increasing 
temperature the phase coherence is lost due to the long- 
wavelength phase fluctuations and the gas may undergo 
an entropy driven, first-order phase transition to the 
Berezinskii-Kosterlitz-Thouless-like phase. This phase 
consists of bound pairs of vortices with opposite circula- 
tion. The gas exhibits superfluidity on both sides of this 
pair creation transition, although above it, only a topo- 
logical order exists in terms of correlations between the 
locations of the vortices. These correlations yield alge- 
braic long-range order, allowing the system to remain in a 
superfluid state. At higher temperatures the paired vor- 
tices unbind at the BKT superfluid to normal phase tran- 
sition. Our present formalism does not allow us to probe 
this upper transition in great detail. Nevertheless, the 
binding energy of a pair is roughly equal to the pair cre- 
ation energy, as is verified by our numerical simulations. 
Therefore, the temperature at which single isolated vor- 
tices can form spontaneously would be twice the temper- 
ature at which bound vortex pairs are created, suggesting 
the temperature range for the BKT phase to be roughly 
equal to that of the phase coherent BEG phase. This esti- 
mation is also in reasonable agreement with the ideal gas 
prediction for the superfluid-normal transition tempera- 
ture Tq. At any rate, the width of the BKT phase should 
exceed the temperature resolution in currently feasible 
experiments. In this context it is worth re-emphasizing 
that in a homogeneous system, BEG only occurs at T = 
and therefore the BKT phase is the only superfluid phase 
of the system at finite temperature. However, in the 
presence of an external confining potential, phase fluctu- 
ations are suppressed, allowing the formation of a phase 
coherent BEG. This leads to the existence of two distinct 
superfluid phases. 

Discussion 

We have calculated the free energy for a dilute har- 
monically trapped two-dimensional Bose gas at ultra-low 
temperatures. We have shown that in such a system 
there may exist two different superfluid phases — a phase 
coherent Bose-Einstein condensate and a phase fluctu- 
ating Berezinskii-Kosterlitz-Thouless-like phase. We ob- 
tain an estimate for the critical temperature separating 
these two phases in terms of the external system parame- 
ters. Since we are interested in the onset temperature for 
spontaneous pair formation, it suffices to only consider a 
single pair of vortices instead of the full interacting pair 
vortex plasma. After this an avalanche of pair creation 
may occur. Unbinding of these vortex pairs eventually 
leads to the destruction of superfluidity in the system. 

The approach adopted here could be extended to take 
the thermal quasiparticle effects explicitly into account, 
but these effects are of higher order and are not expected 
to alter the qualitative picture obtained. This is espe- 



cially true since we are interested in the free energy dif- 
ference rather than its absolute value. In fact, the most 
significant effect neglected in this respect is probably the 
finite-temperature core broadening of vortices due to the 
thermal cloud pressure, which is only significant at high 
temperatures with very small condensate fraction. The 
transition to the BKT phase is predicted to occur at much 
lower temperatures, however, validating our approach. 
The actual transition temperature may also be lower than 
our prediction, since we consider vortex pair creation in 
the vicinity of the trap centre where the self-energy of the 
pair is the highest. In practice one expects the sponta- 
neous pair creation to begin near the low density bound- 
aries of the condensate, which also opens the possibility 
of the co-existence of the BEG and BKT phases. 

The concept of the quasi-condensate is inherent within 
the BKT phase, with the vortex pairs located at the 
boundaries of the locally phase coherent domains. It is 
possible that, at lower temperatures, phase fluctuations 
may already destroy the global phase coherence, yielding 
a quasi-condensate prior to the formation of the BKT 
phase. Thus the pair-creation transition implies a quasi- 
condensate but not necessarily the converse. 

We suggest that it would be possible to directly ob- 
serve the paired vortices in the BKT-like phase. The 
BKT phase may pose observational difficulties even if it 
were present in the system. However, one could possibly 
observe the paired vortices using standard imaging pro- 
cedures after a ballistic expansion of the condensate. An- 
other, theoretically tempting possibility, could be to sep- 
arate the vortices of opposite circulation by breaking the 
time-reversal symmetry. This could be done by rotating 
the condensate below the rotationally activated vortex 
nucleation frequency. Subsequently one could measure 
the residual angular momentum of the system, whose fi- 
nite value would be a clear indication of the pre-existence 
of vortex pairs in the BKT phase, since the direct forma- 
tion of single vortices is forbidden. 
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